Electron-phonon interaction and antiferromagnetic correlations 
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We study effects of the Coulomb repulsion on the electron-phonon interaction (EPI) in a model 
of cuprates at zero and finite doping. We find that antiferromagnetic correlations strongly enhance 
EPI effects on the electron Green's function with respect to the paramagnetic correlated system, 
but the net effect of the Coulomb interaction is a moderate suppression of the EPF Doping leads 
to additional suppression, due to reduced antiferromagnetic correlations. In contrast, the Coulomb 
interaction strongly suppresses EPI effects on phonons, but the suppression weakens with doping. 
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There are strong experimental indications that the 
electron-phonon interaction (EPI) plays a substantial 
role for properties of high-T c cuprates 0, , and that it 
even can lead to formation of small polarons for undopcd 
cuprates 0. The Coulomb interaction in the copper- 
oxide plane is expected to strongly suppress charge fluc- 
tuations. This is often described in the Hubbard or t-J 
models, for which important phonons couple to such 
charge fluctuations ||. One might therefore actually ex- 
pect the Coulomb interaction to strongly suppress the 
EPI. This was indeed found with dynamical mean-field 
theory jdMlcalculations in the paramagnetic phase (P- 
DMFT) yiHlGIl- On the other hand, calculations for 
the undoped antiferromagnetic t-J model using the self- 
consistent Born approximation (SCBA) [llj . or approx- 
imations going beyond the SCBA [l^, found that the 
Coulomb interaction enhances EPI effects on the elec- 
tron Green's function. 

The half-filled Hubbard model becomes an insulator 
for large U. In the P-DMFT the only mechanism for this 
is the quasiparticle weight Z going to zero 0. Such a 
reduction of Z tends to strongly suppress the EPI flllllH , 
which is an important reason for the small effects of the 
EPI in the P-DMFT. By allowing for antiferromagnetism 
(AF) in the DMFT Hi QUI], it is possible to have an 
insulating state although Z remains finite. This suggests 
that it is important to allow for AF when describing the 
EPI, i.e., using an AF-DMFT. This is also suggested by 
the fact that in the SCBA, the enhancement of the EPI 
grows with the value of the exchange constant J ^fjj i.e., 
with the importance of the AF. 

Here we therefore use an AF-DMFT formalism and 
we first consider a half-filled system. In contrast to pre- 
vious work, we find that within the Holstein-Hubbard 
model the effects of the EPI on the electron Green's func- 
tion are neither strongly suppressed nor enhanced by the 
Coulomb interaction. While previous work could only ad- 
dress ^2 the half-filled case, the present formalism 



makes it possible to treat doped metallic cuprates, which 
are of particular interest. We show that the EPI remains 
important for the electron Green's function as long as AF 
is important. Due to the weakening of AF correlations 
as the system is doped, we find that doping reduces the 
effects of the EPI on the electron Green's function and it 
weakens the tendency to polaron formation, in agreement 
with experimental results 

We also calculate the renormalization of the phonon 
frequency. For the undoped system, we find that the 
Coulomb repulsion very strongly suppresses the renor- 
malization. As the system is doped, however, the renor- 
malization of the phonon frequency increases. This is the 
opposite behavior to what we find for the electron Green's 
function. The width and softening of the phonon spec- 
tral function is often used to estimate the strength of the 
EPI ^(| . These results show that for strongly correlated 
systems this approach may strongly underestimate the 

epi (HEi. 

We study the Holstein-Hubbard model 
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where c ia creates an electron with spin a on site % and 
bi a phonon on site i, t > is a hopping integral, U 
the Coulomb interaction between two electrons on the 
same site, ujq is the phonon energy, g a coupling con- 
stant and rii measures the number of electrons on site 
i. We assume an infinite-coordination Bcthc lattice with 
the half band width D and the density of states (DOS) 
N(e) = (2/tvD 2 )^/D 2 — e 2 . We define a dimensionless 
coupling constant A = g 2 /(ujqD). For large U the Hub- 
bard model is approximately equivalent to the t-J model. 
For a two-dimensional (2d) lattice, these models are re- 
lated via J/t = D/U. 
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We solve the DMFT equations for the temperature 
T = 0. The associated impurity problem is solved us- 
ing the Lanczos method. The Hubert space is limited by 
only allowing up to N p h phonons, where N p h ~ 30 de- 
pends on the parameters. The energies of and couplings 
to the bath levels are determined from a continued frac- 
tion expansion for the large U half-filled case and 
otherwise by a fit of the cavity Green's function on the 
imaginary axis |2(ij| . We use up to 25 bath levels. 

We here focus on the quasi-particle weight Z , since our 
criterion for polaron formation is Z being exponentially 
small and since Zq, calculated for A = 0, is expected to 
be crucial for the EPI, as discussed above. Fig. shows 
Zq as a function of U for the half-filled Hubbard model 
according to AF-DMFT and P-DMFT and as a function 
of J for the (2d) t-J model according to the SCBA. The 
SCBA results agree well with exact diagonalization re- 
sults for small clusters [2l|. The half-filled system is an 
insulator in the AF-DMFT and for U > 3D in the P- 
DMFT 0]. While Z drops to zero very quickly with 
U in the P-DMFT, it remains finite in the AF-DMFT. 
For U/D values where the Hubbard and t-J models are 
approximately equivalent, AF-DMFT and SCBA agree 
rather well. This good description of Z suggests that 
AF-DMFT may describe the EPI well. 
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FIG. 1: Quasiparticle weight Zq for A = for the Hubbard 
model (as a function of U/D) according to P-DMFT and AF- 
DMFT and for the 2d t-J model (as a function of J/t) accord- 
ing to the SCBA. The figure illustrates that the AF-DMFT 
gives reasonable values for Zq. 

We first discuss the results for a weak EPI. For nonin- 
teracting electrons (U = 0), the reduction of Z by the 
EPI is given by (l/Z - 1)/A = 4/V for small A and 
uiq <C D, as shown by the arrow in Fig. [5] To determine 
the effect of the EPI for interacting electrons (U > 0), we 
calculate (Zo/Z — 1)/A, shown in Fig. [3 This quantity 
measures how efficiently the EPI reduces the quasiparti- 
cle weight Z with respect to Zq obtained in the absence 
of EPI. In a P-DMFT calculation it was found that the 



EPI very quickly becomes inefficient when U is increased 
[Toj . In contrast, allowing for AF, we find that the EPI 
remains much stronger as U is increased, although it is 
still reduced compared to the noninteracting case. These 
results show that AF is crucial for the EPI of the half- 
filled system. For U/D values where the Hubbard and t-J 
models can be compared, the AF-DMFT and the SCBA 
agree well for loq = 0.025Z?. For larger phonon frequen- 
cies, however, we find that the EPI is appreciably more 
efficient in the SCBA than in the AF-DMFT. 
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FIG. 2: (Zo/Z - 1)/A in the limit A -> for the Hubbard 
model according to the P-DMFT and AF-DMFT and for the 
t-J model according to the SCBA. This quantity indicates 
the effectiveness of the EPI in the weak-coupling limit. The 
figure illustrates how the EPI is much more efficient in the AF- 
DMFT than the P-DMFT, and that AF-DMFT and SCBA 
agree rather well. 

We next focus on strong EPI. Fig. [3] shows Z as a 
function of A for different U . The result for U = 3.5D 
(J/t = 0.29) can be compared with a calculation for the 
2d t-J model 0] (J/t = 0.3 and the same ui ). As A is 
increased, Z is strongly reduced, signaling polaron forma- 
tion. This happens at a somewhat larger critical value A c 
than was found for the t-J model 0], indicated by an ar- 
row. The deviation from Ref. [l2l| is probably mainly due 
to our use of the AF-DMFT and the neglect of "crossing" 
diagrams in Ref. ^3 ■ Good agreement is also found with 
results for the infinite dimension t-J model These 
comparisons suggest that the AF-DMFT is rather accu- 
rate for the half-filled Holstein-Hubbard model. 

For U = 0, Z drops very quickly as a function of A and 
(bi-)polarons are formed at A c w 0.33. As U is increased 
the drop is slightly less rapid and polaron formation hap- 
pens at somewhat larger A c . The Coulomb interaction 
therefore moderately suppresses polaron formation. 

In P-DMFT calculations it was found that the effective 
mass to* depends only weakly on ujq [loj]. In AF-DMFT 
we find a sizable isotope effect on m* and for A ~ A c the 
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FIG. 3: Z as a function of A for different U and for coo — 
0.025D. The arrow shows A c of the t-J model for J/t — 
0.3 (U/D = 3.3) Q. The figure shows how the Coulomb 
interaction moderately suppresses polaron formation [Z — > 
0). 



FIG. 4: Z as a function of A for different fillings n and asso- 
ciated magnetic moments m for U — 3.5D and ujo — 0.025D. 
The figure illustrates how the critical A c is increased as the 
filling is reduced (doping is increased) due to a reduction of 
the the antiferromagnetic correlations. 



effect is comparable to the experimental value |23j . 

Above we have discussed the half-filled system exten- 
sively, since we can compare with other methods and test 
the reliability of the AF-DMFT. The doped cuprates, 
however, are more interesting and challenging, and we 
now focus on them. Fig. 0] shows Z as a function of A 
for U = 3.5D and for different dopings. As the filling 
is reduced (hole doping increased) the staggered magne- 
tization to is reduced. The figure shows how this leads 
to an increase in A c . In a P-DMFT calculation 011 
the other hand, a reduction of the filling leads to a re- 
duction of A c . The increase of A c in the AF-DMFT with 
increased doping is therefore indeed due to the reduction 
of to, since at constant m = 0, A c decreases with doping. 

We find the AF-P transition for U = 3.5D at n = 0.84, 
corresponding to a much larger doping (0.16) than found 
experimentally. This is only partly due to our neglect of 
second nearest neighbor hopping which would introduce 
magnetic frustration in the system. The main reason is 
that in the P state there are also AF correlations which 
lower the energy, but which arc neglected in a DMFT 
calculation. The AF-DMFT calculation therefore favors 
the AF state. To obtained a balanced treatment it is 
necessary to use a cluster DMFT method 0, . Such 
a calculation would introduce AF correlations also in the 
P state, and like in the AF-DMFT calculation these cor- 
relations would weaken as the doping is increased. This 
should increase A c with doping in a qualitative similar 
way as in Fig. 

Experimentally, polaron formation is found to disap- 
pear as the system is moderately doped Q . This may be 
partly due to screening of the EPI, leading to a reduction 
of A. However, the suppression of polaron formation with 



doping for fixed A, illustrated in Fig.0] should also be an 
essential part of the explanation. 

To study effects of the Coulomb interaction, earlier 
work compared with a Holstein model with a single elec- 
tron at the bottom of the band [i"llll2j |. A better compar- 
ison, however, is with a half-filled Holstein model, since 
we can then increase U keeping the number of electrons 
unchanged. For a half-filled 2d Holstein model such a 
comparison is difficult, due to Fermi surface nesting and 
phase transitions. In DMFT calculations using a Bcthe 
lattice this is not a problem. 

The EPI appears much stronger for the half-filled Hol- 
stein model than for a single electron at the bottom of the 
band. Comparing with the half-filled case, we therefore 
find that the EPI is suppressed by the Coulomb inter- 
action, while earlier studies, comparing with the single 
electron case, found an enhancement. 

To discuss the difference between the two cases for a 
weak EPI, we calculate the electron-phonon part of the 
electron self-energy (g and u> small, uiq <C D) 

2a 2 

RcI] cp H = -a— y — u, (2) 

where a = 1 for the 2d single electron case but a = 4 
for the half-filled case with a semi-elliptical DOS. This 
large difference is partly due to the DOS of the 2d Hol- 
stein model being smallest at the bottom of the band, 
and partly due to Re E op having contributions both from 
higher and lower states for the half-filled case, but only 
from higher states for the single electron case. 

To understand the difference for a strong EPI, we study 
polaron formation in the adiabatic limit by comparing 
states with free electrons and perfectly localized electrons 
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[26j. We find Ef ree = —Apt per electron, where /? = 1 
(one electron) or (3 — 4/(37r) ss 0.42 (half-filled case), 
and E\ oc = — g 2 / w ph per electron for both cases. We 
assume that polarons form when \E\ oc \ > |i?f rco |. This 
leads to a large A c = 1 [53 for a single electron and a 
much smaller A c = 0.42 for the half-filled case |29j . 

We emphasize the remarkably small value, A c = 0.33, 
for polaron formation in a half-filled Holstein model, 
meaning that Migdal's theorem breaks down for quite 
small A. Using different values for ujq, small A c have also 
been obtained earlier jMUl HIIH using DMFT calcu- 
lations. 
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FIG. 5: Phonon spectral function for different values of A. 
The bare phonon frequency is uoo = 0.LD and a Lorentzian 
broadening with the full width half maximum of 0.04 D has 
been introduced. The figures a-c show how the phonon soft- 
ening at half-filling is dramatically suppressed by U and figure 
d that the softening increase with doping S. 

We now consider the influence of the EPI on phonons. 
Fig.|SJi shows the phonon spectral function for a small U 
at half-filling (doping (5 = 0) for different EPI strengths. 
The figure illustrates how the phonon is softened sub- 
stantially as the EPI is increased. Fig. [5j>c show this 
softening is strongly reduced when U is increased. The 
reason is that for large U charge fluctuations are strongly 
suppressed, and the system can only respond weakly to 
a phonon which couples to the net charge on the atoms. 
Fig. [SJi shows how the softening increases as the dop- 
ing is increased, due to the doped holes responding to 
phonons. This is in agreement with neutron scattering 
measurements 0. The figure illustrates that the influ- 
ence of the EPI on the phonon self-energy is dramatically 
different from the influence on the electron self-energy. 
The reason is that the electron self-energy measures the 
response of the system to the removal or addition of a 
charge, which leads to a strong response even in cases 
where charge fluctuations are otherwise suppressed. 

While paramagnetic DMFT calculations for the 



Holstcin-Hubbard model show that effects of the EPI 
on electrons (quasiparticle weight) arc very strongly sup- 
pressed by the Coulomb interaction, we find that this 
suppression is only moderate when antiferromagnetic 
(AF) correlations are included. As the doping is in- 
creased, the AF is reduced and the EPI is more sup- 
pressed. In contrast, at half-filling, the Coulomb inter- 
action strongly suppresses effects of the EPI on phonons 
(phonon softening) , while the suppression is reduced with 
doping. These trends are consistent with experiment. 
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